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We study the time variation of fundamental constants in the early Universe. Using data from 
primordial light nuclei abundances, CMB and the 2dFGRS power spectrum, we put constraints 
on the time variation of the fine structure constant a, and the Higgs vacuum expectation value 
< u > without assuming any theoretical framework. A variation in < t; > leads to a variation in 
the electron mass, among other effects. Along the same line, we study the variation of a and the 
electron mass m^. In a purely phenomenological fashion, we derive a relationship between both 
variations. 
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C/3 , I. INTRODUCTION 

Unification theories, such as super-string [l|, 0, H, 0, H, 0], brane- world 0, H, H, and Kaluza-Klein theories 
[ill [l^ . [isl . O ITsl l. allow fundamental constants, such as the fine structure constant a and the Higgs vacuum 
expectation value < t; >, to vary over cosmological timescales. A variation in < t; > leads to a variation in the electron 
mass, among other eff'ects. On the other hand, theoretical frameworks based in first principles, were developed by 
different authors [l6l . [itI . [18| in order to study the variation of certain fundamental constants. Since each theory 
1 predicts a specific time behaviour, by setting limits on the time variation of fundamental constants some of these 
theories could be set aside. 

Limits on the present rate of variation of a and ji = ^ (where rrie is the electron mass and rup the proton mass) 

00 ! are provided by atomics clocks [l^ HH, [H, [H, [23| . Data from the Oklo natural fission reactor [2^, [2^ and half lives 
■ of long lived (3 decayers [13] allow to constrain the variation of fundamental constants at z ~ 1 . Recent astronomical 
data based on the analysis of spectra from high-redshift quasar absorption systems suggest a possible variation of a 
and ^ [2i, [1^, [13, [31, 32, 33, 34]. However, another analysis of similar data gives null variation of a [35, 36, 37, 38]. 
Big Bang Nucleosynthesis (BBN) and Cosmic Microwave Background (CMB) also provide constraints on the variation 
of fundamental constants. Although the limits imposed by BBN and CMB, are less stringent than the previous ones, 
they are still important since they refer to the earliest cosmological times. 

In previous works, we have studied the time variation of the fine structure constant in the early Universe to 
test Bekenstein model t39j] and the time variation of the electron mass to test Barrow-Magueijo model 40]. However, 
unifying theories predict relationships among the variation of gauge coupling constants which depend on the theoretical 
framework. In this work, we perform a phenomenological analysis of the joint time variation of a and < u > in the 
early Universe without assuming a theoretical framework. 

The model developed by Barrow & Magueijo [l^ predicts the variation of rrie over cosmological timescales. This 
model could be regarded as the low energy limit of a more sophisticated unified theory. In such case, the unifying 
theory would also predict variation of gauge coupling constants and in consequence the variation of a. Thus, in order 
to provide bounds to test such kind of models, we also analyze in this paper the joint variation of a and mg without 
assuming a theoretical framework. 

The dependence of the primordial abundances on a has been analyzed by Bergstrom et al. 14111 and improved by 
NoUet & Lopez 4?], while the dependence on < > has been analyzed by Yoo and Scherrer [43l]. Semi-analytical 
analyses have been performed by some of us in earlier works [3, [4^ . Several authors [i^, [13, l48l] studied the effects 
of the variation of fundamental constants on BBN in the context of a dilaton superstring model. Miiuller et al 
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pg | calculated the primordial abundances as a function of the Planck mass, fine structure constant, Higgs vacuum 
expectation value, electron mass, nucleon decay time, deuterium binding energy, and neutron-proton mass difference 
and studied the dependence of the last three quantities as functions of the fundamental coupling and masses. Coc et 
al. ^50*1 set constraints on the variation in the neutron lifetime and neutron-proton mass difference using the primordial 
abundance of '*He. Cyburt et al. 51] studied the number of relativistic species at the time of BBN and the variations 
in fundamental constants a and Gn- Dent et al [52| studied the dependence of the primordial abundances with 
nuclear physics parameters such as Gjv, nucleon decay time, a, rrie, the average nucleon mass, the neutron-proton 
mass difference and binding energies. Finally, limits on cosmological variations of a, A^qcd and quark mass (to*) 
from optical quasar absorption spectra, laboratory atomic clocks and from BBN have been established by Flambaum 
et al. 

In this paper, we study the effects of a possible variation of a and < f > on the primordial abundances, including 
the dependence of the masses of the light elements on the cross sections, and using the dependence on < w > of the 
deuterium binding energy calculated in the context of usual quantum theory with a phenomenological potential. We 
use all available observational data of D,'''He and ^Li to set constraints on the joint variation of fundamental constants 
at the time of BBN. 

However, we do not consider a possible variation of A^qcd- Indeed, the dependence of the physical quantities 
involved in the calculation of the primordial abundances with a varying Aqcd is highly dependent on the model. 
The analyses of refs. [i^, [13], for example, are done in the context of a string dilaton model. Therefore, we will 
not consider such dependencies even though it has been analyzed in the literature [46, 47, 53, 54, 55]. Our analysis, 
instead, is a model independent one. 

Previous analysis of CMB data (earlier than the WMAP three- year release) including a possible variation of a have 
been performed by refs. 5^ 53, Is^ and including a possible variation of rUe have been performed by refs. (43l . IsF ^ 



The work of Ichikawa et al. |58j is the only one that assumes that both variations are related in the context of string 
dilaton models. In this work, we follow a completely different approach, by assuming that the fundamental constants 
vary independently. 

The paper is organized as follows. In section |TT1 we present bounds on the variation of the fine structure constant 
and the Higgs vacuum expectation value during Big Bang Nucleosynthesis. We also discuss the difference between 
considering < v > variation and TOg variation during this epoch. In section Hill we use data from the CMB and from 
the 2dFGRS power spectrum to put bounds on the variation a and < v > (or nie) during recombination, allowing also 
other cosmological parameters to vary. In section HVl we use the a — nif. and a— < v > confidence contours to obtain 
a phenomenological relationship between both variations and then discuss our results. Conclusions are presented in 
section IVl 



II. BOUNDS FROM BBN 



Big Bang Nucleosynthesis (BBN) is one of the most important tools to study the early universe. The standard model 
has a single free parameter, the baryon to photon ratio r]B, which can be determined by comparison between theoretical 
calculations and observations of the abundances of light elements. Independently, the value of the baryonic density 
Q.sh?' (related to tjb) can be obtained with great accuracy from the analysis of the Cosmic Microwave Background 
data [53, |6(i, i61j . Provided this value, the theoretical abundances are highly consistent with the observed D but not 
with all ^He and ^Li data. If the fundamental constants vary with time, this discrepancy might be solved and we may 
have insight into new physics beyond the minimal BBN model. 

In this section, we use available data of D, "^He and ^Li to put bounds on the joint variation of a and < v > and on 
the joint variation of a and rrie at the time of primordial nucleosynthesis. The observational data for D have been taken 
from refs. [H, M, M, M, M, M, SJ^LFor ^Li we consider the data reported by refs. [zO, [Zll, Izl, III, III Izl Izl ■ 
For ^He, we use the data from refs. [TTllTSj (see ref.[3§] for details). 

We checked the consistence of the each group of data following ref. [t^ and found that the ideogram method 
plots are not Gaussian-like, suggesting the existence of unmodelled systematic errors. We take them into account by 
increasing the errors by a fixed factor, 2.10, 1.40 and 1.90 for D, *He and ''Li, respectively. A scaling of errors was 
also suggested by ref. [80] . 

The main effects of the variation of the fine structure constant during BBN are the variation of the neutron to 
proton ratio in thermal equilibrium produced by a variation in the neutron-proton mass difference, the weak decay 
rates and the cross sections of the reactions involved during the first three minutes of the Universe. The main effects 
of the variation of the Higgs vacuum expectation value during BBN are the variation of the electron mass, the Fermi 
constant, the neutron-proton mass difference and the deuterium binding energy, affecting mostly the neutron to proton 
ratio, the weak decay rates and the initial abundance of deuterium. In appendix [X] we give more details about how 
the physics at BBN is modified by a possible change in a, < w > and TOg. We modify the Kawano code [8l[ in order 
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to consider time variation of a and < v > and time variation of a and rUe during BBN. The coulornb, radiative and 
finite temperature corrections were included following ref. (s^l • We follow the analysis of refs. [4l|, \^ to introduce 
the variation in a on the reaction rates. The main effects of a change in a in nuclear reaction rates are variations in 
the Coulomb barrier for charged-induced reactions and radiative captures. We introduce the dependence of the light 
nuclei masses on a, correction that affect the reaction rates, their inverse coefficients and their Q- values [ij]. We also 
update the value of the reaction rates following ref. [4l| . 

To illustrate the effect of the variation in the fine structure constant on the reactions rates, we present the nuclear 
reaction rate ofd + d^p + t (^R[dd;pt] — 0.9310^'^f2B/i^rg TVa < av >) as a function of (Aa = a — ao and ao 
is the current value of the fine structure constant): 



R[dd;pt] 



2.36910-3nB/i^Tj/^ao^/^ ( 1 



Aa 
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Tq'^ - 3.96610V^^^ao^/^ 1 + 



(1) 



where Tg is the temperature in units of 10^ K and /i is the reduced mass. The reduced mass also changes if the fine 
structure constant varies with time (see appendix]^. This nuclear reaction is important for calculating the final 
deuterium abundance since this reaction destroys deuterium and produces tritium which is crucial to form ^He. In 
Figure [T] we present the value of Na < av > for this reaction as a function of the temperature, for different values of 
If the fine structure constant is greater than its present value, the reaction rate is lower than in the case of no a 
variation. A decrease in the value of this reaction rate results in an increase in the deuterium abundance. 




4 6 
T[10^ K] 



FIG. 1: Na < <jv > (in units of lO'^ cm'' s ^ g ^) for the reaction d-|-d^p-|-t, as a function of the temperature (in units of 
10'' K), when ^ = -0.1 (solid line), ^ = 0.0 (dashed line) and ^ = 0.1 (dotted line) 

The "^He abundance is less sensitive to changes in the nuclear reaction rates than the other abundances (deuterium 
and ^Li) 41] and very sensitive to variations in the parameters that fixed the neutron-to-proton ratio. In thermal 
equilibrium, this ratio is: 

^ = e--^"-/^ , (2) 
Yp 

where Yn (Yp) is the neutron (proton) abundance, Am„p is the neutron-proton mass difference and T is the temperature 
in MeV. When the weak interaction rates become slower than the Universe expansion rate the neutron-to-proton ratio 
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freezes-out at temperature Tj. Afterwards, nearly all the available neutrons are captured in "^He [4l|, this abundance 
can be estimated by: 



I4 



(3) 



where ^ — e '^'""p/^/ . The neutron-proton mass difference is affected by a change in the fine structure constant 
and in the Higgs vacuum expectation value: 

^ = -0.587^ + 1.587^^. (4) 

An increase in the fine structure constant results in a decrease in Ar7i„p, this produces a larger equilibrium neutron- 
to-proton ratio and a larger abundance of ^He. However, an increase in < t; > leads to an increase in Aninp- This 
produces a smaller equilibrium neutron-to-proton equilibrium ratio and a smaller abundance of ^He [43| . 

The freeze-out temperature of weak interactions is crucial to determinate the amount of available neutrons and 
therefore the primordial abundance of ^He. This temperature is modified if the Higgs vacuum expectation value is 
changed during BBN due to changes in the weak reaction rates (see appendix!^. A larger Higgs vacuum expectation 
value during BBN results in: i) a smaller Gp leading to earlier freeze-out of the weak reactions (n p), producing 
more ^He; ii) an increase in We, a decreasing of n <-> p reaction rates and also producing more "^He (isj . 

The dependence of the deuterium binding energy on the Higgs vacuum expectation value is extremely model 
dependent. Bean and Savage [s^l studied this dependence using chiral perturbation theory and their results were 
applied by several authors |43i . i49l. |55il . We performed another estimation in the context of usual quantum theory, 
using the effective Reid potential [8J] for the nucleon-nucleon interaction (paper in preparation). Even though Yoo 
and Schcrrer had shown that very different values for lead to similar constraints on the change of < v >, we 
perform our calculation using two different relationships (see Table |l|: i) the obtained by ref. [1^ using the results 
of ref. [13] ; ii) the obtained using the effective Reid potential. From Table HI it follows that the value obtained using 
the effective Reid potential lies in the range allowed by the estimation of Beane and Savage [s^ . The variation of the 
deuterium binding energy due to a time variation of the Higgs vacuum expectation value is related to as: 

Aeji _ dep A <v > 

{€d)q 9m^2(eD)o < v >o 

We call K = ,7"', hereafter. 



TABLE I: Values used in this work for and the coefficient k in the relationship f^'? = ^ ^<"> , 





dep 
dm„ 


K 


Yoo and Scherrer 


-0.159 


-5.000 


Reid potential 


-0.198 


-6.230 



An increase in the Higgs vacuum expectation value results in a decrease in the deuterium binding energy, leading 
to an smaller initial deuterium abundance: 

Y V „11.605eD/T9 

Yd = TT^ , (6) 

0.47110-iorg' 

where ep is in MeV. The production of ^He begins later, leading to a smaller helium abundance but also to an increase 
in the final deuterium abundance [43 |. 

To assume time variation of the electron mass during BBN is not exactly the same as to assume time variation of 
the Higgs vacuum expectation value since the weak interactions are important during this epoch. There exist some 
well tested theoretical model that predict time variation of the electron mass [l^ . For this reason we compute the 
light nuclei abundances and perform a statistical analysis using the observational data mentioned above to obtain the 
best fit values for the parameters for the following cases: 



• variation of a and < v > allowing rjs to vary. 



5 



• variation of a and < v > keeping r]B fixed, 

• variation of a and rUe allowing rjB to vary, 

• variation of a and me keeping r]B fixed. 

Even though the WMAP data are able to constrain the baryon density with great accuracy, there is still some 
degeneracy between the parameters involved in the statistical analysis. For this reason, we allow the joint variation 
of baryon density and the other two constants to obtain an independent estimation for tjb- In the cases were tjb is 
fixed, we assume the value reported by the WMAP team ^ {vb^'^^^ = (6.108 ± 0.219)10-^°). We also present in 
Table |lTl as an independent estimation, the best fit for the baryon density when all the constants are fixed at their 
present value. It is shown that the only reasonable fit is found when the lithium data is removed from the data set, 
and the value for rjB is consistent with the value reported by Spergel et al. (60|. The fits considering the data of ^Li 
are not reasonable and the best fit for rjB is not consistent with the WMAP value. 

TABLE II: Best fit parameter values, la errors for the BBN constraints on rys (in units of 10^^") keeping all the fundamental 
constants fixed at their present value. 





77s ± a [10""'" 


JV-1 


D He +^ Li 


4.310 ±0.050 


10.00 


*He +^ Li 


3.920 ± 0.080 


6.53 


D+'^Li 


4.270 ± 0.060 


9.27 


D+* He 


o- ' ru_Q 


1.61 



Table Hill shows the results for the analysis of the variation of a and < v > when 773 is allowed to vary. These results 
correspond to the relationship between cd and < v > obtained using Reid potential. This fit is consistent within 
1(7 with the acquired considering the value of k calculated by Yoo and Scherrer. We will not perform the statistical 
analysis again excluding one group of data for the case where 775 is allowed to vary since we would have two groups 
of data and three unknown variables. A reasonable fit can be found when a, < w > and rjB are allowed to vary. This 
fit is consistent within 6f7 with non null values for the variations of a and < v > while the obtained value for i]b is 
not consistent with the estimation of WMAP within 3(7. 

Left Figure [2] shows a strong degeneracy between and 775, ^^Z^ and rjB and ^ and ^^Z^ ■, (A < u >=< 
u > — < 1; >o and < u >o is the present value of the Higgs vacuum expectation value). From this phenomenological 
approach, the variation of the fundamental constants can be used to reconcile the observed primordial abundances. 

TABLE HI: Best fit parameter values and Icr errors for the BBN constraints on — , ^<"> /^™<= allowine riB fin units of 

^ CVQ ' <v>Q ' (me)Q ' o /-D V 

10""') to vary and considering D He +^ Li. We use the estimation 7^^^^- = -6.230#^^ to obtain the fit on and 

1 J fall (Er>)o <v>o CO ' <">o 

riB- 





^<^> ± a 

<v>o 




r7B±a[10-^»] 


JV-3 


0.198ini^ 






8.005l°;^g 


1.16 


0.210lHi^ 




-0 250+''''"5 


7 rqq+0.447 

' •>J'J'J-0.5O2 


1.11 



When a, rrie and rjB are considered as free parameters, we also obtain a reasonable fit (see Table HIT)) . Once again, 
we will not perform the statistical analysis again excluding one group of data since we would have two groups of data 
and three unknown variables. There is consistency within 6(7 with variation of a and mg but the obtained value for 
rjB is not consistent with the estimation of WMAP within 3(7. Right Figure [2] shows a strong degeneracy between 
and rjB, and rjB and ^ and (Amg = — (?7ie)g, {me)g is the present value of the electron mass). 

From this result it might be possible to reconcile the observed primordial abundances and the WMAP estimation for 
the baryon density. 

By comparing the results presented in Table IIIII and in Table |TT1 it can be noticed that the variation of the 
fundamental constants improves the statistical analysis providing a Xmin / ~ value that is closer to one. 

We now consider the joint variation of the fine structure constant and the Higgs vacuum expectation value with rjB 
fixed at the WMAP estimation. In this case, it is reasonable to repeat the analysis excluding one group of data at the 
time. The results are presented in Table ITVl and were obtained using ^y- = — 6.230Ag^. The fits obtained using 
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FIG. 2: Left Figure: Likelihood contours for ^Kvy^ ^'i'^ '7s (in units of 10 ^'^) and 1 dimensional Likelihood (using 

= -6.230^^). Right Figure: Likelihood contours for and ris (in units of 10"^°) and 1 dimensional 

Likelihood. 



= —5.000-=^^ are consistent, within Itr, with the ones presented. There is good fit for the whole data set and 
also excluding one group of data at each time. In any case, there is a strong degeneracy between ^ and (see 
Figure Considering all data or *He +^ Li we find variation of both fundamental constants, a and < t; >, even at 
6(7. However, if the statistical analysis is performed with D +^ He we find null variation for both constants within Icr. 
For completeness we also modified the Kawano's code in order to calculate the different primordial abundances for 

two values inside the range of calculated by refs. [13, [s^ ^—0.15 < < —0.05^ and performed the statistical 

test in order to obtain the constraints on the variation of the fundamental constants a and < w >. All the results are 
consistent within 2(7 with the ones presented above. 



TABLE IV: Best fit parameter values, la errors for the BBN constraints on ^ and with rjB fixed at the WMAP 

estimation. The results correspond to the estimation , '^'P = —6.230-%^^. 



Data 


an 


^<:> ± a 
<v>n 


N-2 


D He Li 


0.140 ±0.006 


0.032 ± 0.002 


2.52 


■^He +^ Li 


n i/i«+''-004 

U-f48_o.008 




U.UOO_Q 


1.23 


Li 


090+" °^^ 


-0 070+° °^" 

U.U/U_o 026 


1.15 


D+* He 


U.UOU_Q Q3Q 


_0 002+° "°'' 

U.UUZ_Q 008 


1.03 



Table |V] and Figure [4] show the results obtained when only a and mg are allowed to vary. There is a strong 
degeneracy between the variations of a and the variations of me in all the cases considered. We find reasonable fits 
for the whole data set and also excluding one group of data at each time. Considering all data or ^He +^ Li, we find 
variation of a and me, even at 6a. On the other hand, if the statistical analysis is performed with D +^ He we find 
null variation for both constants within la. 

Richard et al. [s^ have pointed out that a better understanding of turbulent transport in the radiative zones of 
the stars is needed in order to get a reliable estimation of the ^Li abundance, while Menendez and Ramirez \sf\ have 
reanalyzed the ^Li data and obtained results that are marginally consistent with the WMAP estimate. On the other 
hand, Prodanovic and Fields [s^ put forward that the discrepancy with the WMAP data can worsen if contamination 
with ^Li is considered. Therefore, we adopt the conservative criterion that the bounds on the variation of fundamental 
constants obtained in this paper are those where only the data of D and ^He are fitted to the theoretical predictions 
of the abundances. This paper shows evidence for variation of fundamental constants in the early Universe if the 
reported values for the ^Li abundance are confirmed by future observations and/or improvement of the theoretical 
analyses. 

In Table IVTl we summarize our results for the variation of a, < v > and rUe, using D +^ He data in the statistical 
analysis. The sign of a variation is the same for all the cases. However, the sign of the variation of < f > or mg 
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FIG. 3: Likelihood contours for — vs ^^"^ and 1 dimensional Likelihood. Left Figure: D +^ Li data; Rieht Figure: D +■* He. 



The Figures were obtained using 



(«r>)o 



-6.230^. 



TABLE V: Best fit parameter values, Icr errors for the BBN constraints on and 
estimation. 



(me)o 



with r\B fixed at the WMAP 



Data 




Ame 1 ^ 

(™c)n 


iV-2 


D +^ He +^ Li 


0.159 ±0.008 


-0.213 ±0.012 


1.85 


■'He +^ Li 


0.163 ±0.008 


-0.218 ± 0.013 


1.00 


0+^" Li 


0.067i°;°o?^ 


0.447 ± 0.134 


1.00 


0+" He 


-0.036l°;°g 


020+° °^® 

U.UZU_() Qg4 


1.00 




0.14 0.15 0.16 0.17 0.18 -0.25 -0.23 -0.21 -0.19 -0.17 
Aot/oto Amg/(mg)o 




-0.2 -0.1 
Aa/ao 



0.15 -0.05 0.05 0.15 0.25 
Ame/(me)o 



FIG. 4: Likelihood contours for — vs , '"f and 1 dimensional Likelihood. Left Figure: D He +' Li data; Right Fi gure: 



D +^ He. 



changes depending on whether the joint variation with a is considered or not. In all cases, we have found null variation 

2 

of the fundamental constants within Scr and the values of -^f^ (where g = 2 for the ease where two constants are 
allowed to vary, and g = I when only one fundamental constant is allowed to vary) are closer to one, resulting in 
reasonable fits for all the cases. 
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TABLE VI: Comparison between the best fit parameter values and la errors for the BBN constraints on 



An A<u> 



and 



Ame 



ag ' <V>Q (™c)o ' 

with q — 2 for the case where two 



when one or two constants are allowed to vary [39|, |40|. We also present the values of -j^^ 
constants are allowed to vary, whereas g = 1 when only one fundamental constant is allowed to vary. We consider D +^ He 
data and r^B fixed at the WMAP estimation. 





Time variation of 


Parameter 


Q and < V > 


a and rUe 




<v > 


rrie 


an 


-0 nsn+° °^5 

U.UOU_Q 030 


-0 036+" °^2 


-0.020 ± 0.007 






t<l> ± a 

<v>n 


-0 002+° °°'' 

U.UUZ_o 008 






0.004 ± 0.002 




Ame _L ^ 
(™-e)n 




U.UZU_o.o64 






-0.024 ± 0.008 




1.03 


1.00 


0.90 


0.97 


0.95 



III. BOUNDS FROM CMB 

The cosmological parameters can be estimated by an analysis of the Cosmic Microwave Background (CMB) radi- 
ation, which gives information about the physical conditions in the Universe just before decoupling of matter and 
radiation. 

The variation of fundamental constants affects the physics during recombination (see appendix [B] for details). At 
this stage of the Universe history, the only consequence of the time variation of < w > is a variation in Wg. The main 
effect of a and me variations is the shift of the epoch of recombination to higher z as a or me increases. This is easy 
to understand since the binding energy _B„ scales as a^rrie, so photons should have higher energy to ionize hydrogen 
atoms. In Figs. [5] we show how the ionization history is affected by changes in a and in mg, in a flat universe with 
cosmological parameters [VLbh"^ ,^lcDMh^ ih^r) = (0.0223,0.1047,0.73,0.09). When a and/or me have higher values 
than the present ones, recombination occurs earlier (higher redshifts). The ionization history is more sensitive to a 
than to me because of the B„ dependence on this constants. 



x"" 0.6 




1.2 



0.6 




me/(mJo=1.05 
mg/(mg)o=1 .00 
mg/(mg)o=0.95 



500 



FIG. 5: Ionization history as a function of redshift, for different values of a (left panel) and me (right panel) at recombination 
time. 



The most efficient thermalizing mechanism for the photon gas in the early universe is Thomson scattering on free 
electrons. Therefore, another important effect produced by the variation of fundamental constants, is a shift in the 
Thomson scattering cross section gt, which is proportional to m~^a^. 

The visibility function, which measures the differential probability that a photon last scattered at conformal time 
•q, depends on a and mg. This function is defined as 



. dn 
dr] ' 



where 



dn 
dr] 



(7) 
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FIG. 6: Visibility function as a function of conformal time in Mpc, for different values of a (left panel) and rrie (right panel). 



is the differential optical depth of photons due to Thomson scattering, Up is the total number density of protons (both 
free and bound), Xe is the fraction of free electrons, and a is the scale factor. The strongest effect of variations of a 
and TOe on the visibility function occurs due to the alteration of the ionization history Xe{f])- In Fig. [6] we show that 
if a and/or me were smaller (larger) at recombination than their present values, the peak in the visibility function 
would shift towards smaller (larger) redshifts, and its width would slightly increase (decrease). 

The signatures on the CMB angular power spectrum due to varying fundamental constants are similar to those 
produced by changes in the cosmological parameters, i.e. changes in the relative amplitudes of the Doppler peaks and 
a shift in their positions. Indeed, an increase in a or rue leads to a higher redshift of the last-scattering surface, which 
corresponds to a smaller sound horizon. The position of the first peak (^i) is inversely proportional to the latter, so a 
larger results. Also a larger early integrated Sach- Wolfe effect is produced, making the first Doppler peak higher. 
Moreover, an increment in a or rrie decreases the high-^ diffusion dampin g, w hich is due to the finite thickness of the 
last-scattering surface, and thus, increases the power on very small scales |89l lool. [9l| . All these effects are illustrated 
in Figs. El 
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FIG. 7: The spectrum of CMB fluctuations for different values of a (left panel) and me (right panel). 

To put constraints on the variation of a and < v > during recombination time, we performed a statistical analysis 
using data from the WMAP 3-year temperature and temperature-polarization power spectrum [60l | , and other CMB 
experiments such as CBI [11], ACBAR [ol], and BOOMERANG HIhI], and the power spectrum of the 2dFGRS [11]. 
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We consider a spatially-flat cosmological model with adiabatic density fluctuations, and the following parameters: 

P= (nBh^,ncDMh\e,T,—,^^^,n,,A,] , (8) 
V ao <v>o ) 

where ^cdm^ is the dark matter density in units of the critical density, 6 gives the ratio of the comoving sound 
horizon at decoupling to the angular diameter distance to the surface of last scattering, r is the reionization optical 
depth. Us the scalar spectral index and As is the amplitude of the density fluctuations. 

The parameter space was explored using the Markov Chain Monte Carlo method implemented in the CosmoMC code 
of ref. [92] which uses CAMB ^] to compute the CMB power spectra and RECFAST [9^ to solve the recombination 
equations. We modified these numerical codes in order to include the possible variation of a and < v > (or me) 
at recombination. We ran 8 Markov chains and followed the convergence criterion of ref. jlOCll | to stop them when 
R-\< 0.0149. Results are shown in Table |VlI] and Figure [i 
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FIG. 8: Marginalized posterior distributions obtained with CMB data, including the WMAP 3-year data release plus 2dFGRS 
power spectrum. The diagonal shows the posterior distributions for individual parameters, the other panels shows the 2D 
contours for pairs of parameters, marginalizing over the others. 

It is noticeable the strong degeneracies that exist between ^^^^ and Q^cdm^ ^ ^^"^ and O, — and Uc, and 
also between — and -^f^- The values obtained for ils/i^, /i, VtcDAih^, t, and agree, within la, with those of 
WMAP team |60|, where no variation of a nor < w > is considered. It is interesting to note that our results for the 
cosmological parameters are similar to those obtained considering the variation of one constant at each time [39, 40]. 
Our results are consistent within la with no variation of a and < w > at recombination. 

In Figures [9] and [10] we compare the degeneracies that exist between different cosmological parameters and the 
fundamental constants when one or both constants are allowed to vary. In any case, the allowable region in the 
parameter space is larger when both fundamental constants are allowed to vary. This is to be expected since when the 
parameter space has a higher dimension the uncertainties in the parameters are larger. The correlations of with 
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TABLE VII: Mean values and la errors for the parameters including a and < v > variation. For comparison, results where 
only one fundamental constant is allowed to vary are also shown [39l. |40|]. Hq is in units of km s~^ Mpc~^. 



Parameter 


a and < v > variation 


a variation 


< V > variation 




0.0218 ±0.0010 


0.0216 ± 0.0009 


0.0217 ±0.0010 




0.106 ±0.011 


0.102 ± 0.006 


0.101 ±0.009 


e 


1 

l.UOO_Q 


1.021 ± 0.017 


1.020 ±0.025 


r 


0.090 ± 0.014 


0.092 ± 0.014 




Aa/ao 


-0.023 ±0.025 


-0.015 ± 0.012 




A < i; > / < ?; >o 


0.036 ± 0.078 




-0.029 ± 0.034 


Us 


0.970 ± 0.019 


0.965 ±0.016 


0.960 ±0.015 


As 


3.054 ± 0.073 




3.020 ± 0.064 


Ho 


70Atti 




68.it^:^ 



the other cosmological parameters change sign when < v > is also allowed to vary. On the contrary, the correlations 



of 



A<v> 

<V>0 



with cosmological parameters do not change sign when a is also allowed to vary. 




0.025 




FIG. 9: la and 2a contour levels. Dotted line: variation of a and < v >; solid line: only a variation. The cosmological 
parameters are free to vary in both cases. 
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FIG. 10: la and 2a contour levels. Dotted line: variation of a and < v >; solid line: only < v > variation. The cosmological 
parameters are free to vary in both cases. 

When only one fundamental constant is allowed to vary, the correlation between this constant and any particular 
cosmological parameter has the same sign, no matter whether the fundamental constant is a or < w >. This is because 
both constants enter the same physical quantities. However, since the functional forms of the dependence on a and 
< V > are different, the best fit mean values for the time variations of these fundamental constants are different and 
the probability distribution is more extended in one case than in the other. Nevertheless, in the cases when only one 
constant is allowed to vary, it prefers a lower value than the present one. 
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IV. DISCUSSION 



In section ini we obtained bounds on the variation of a and < v > using the observational abundances of D, ^He 
and ^Li. We performed different analyses: i) we allow ry^ to vary and ii) we keep rjB fixed. We also performed the 
same analyses for two different estimations of the dependence of the deuterium binding energy on the pion mass or 
the Higgs vacuum expectation value: i) the obtained by Yoo and Scherrer who considered the coefficient for the linear 
dependence obtained by Bean and Savage [83]; ii) the obtained using the Reid potential for the description of the 
nucleon-nucleon interaction and without using chiral perturbation theory. The best fits for these two different cases 
are consistent within Icr. Wc found reasonable fit for the variation of a, < v > and rjs for the whole data set and for 
the variation of a, < v >, keeping r]B fixed, for the whole data set and also when we exclude one group of data. We 
only found variation of the fundamental constants when the ''Li abundance is included in the statistical analysis. We 
also calculated the light abundances, keeping rjB fixed at the WMAP estimation, for different values of the dependence 
of on the Higgs vacuum expectation value, inside the range proposed in ref. [85|, and performed the statistical 
analysis. These results are consistent within la with the ones presented in section [Til 

We also considered the joint variation of a and with ijb variable and fixed at the WMAP estimation. In this 
case, we also obtained reasonable fits for the whole data set. When the ^Li abundance was included in the fit, we 
obtained results consistent with variation of fundamental constants (and ijb consistent with the WMAP value). From 
a phenomenological point of view, to vary a and nie solves the discrepancy between the ^Li data, the other abundances 
and the WMAP estimate. However, it is important to mention that the theoretical motivations for me being the 
varying fundamental constant are weak. 

We have discussed in section |TT] that there is still no agreement within the astronomical community in the value 
of the ^Li abundance. We think that more observations of ^Li are needed in order to arrive to stronger conclusions. 
However, if the present values of ''Li abundances are correct, we may have insight into new physics and varying 
fundamental constants would be a good candidate for solving the discrepancy between the light elements abundances 
and the WMAP estimates. 

In section Hill we calculated the time variation of a and < v > (or m^) with data from CMB observations and the 
final 2dFGRS power spectrum. In this analysis, we also allowed other cosmological parameters to vary. We found no 
variation of a and < v > within Icr, and the values for the cosmological parameters agree with those obtained by ref. 
[60| within Icr. 

In Figure [Tl] we show the 2D contours for — and ^^Z'^ obtained from BEN and CMB data. The correlation 

° ' ' ao <v>o 

coefficients are —0.82 for CMB and 0.77 for BBN. There is a small region where the two contours superpose which 
is consistent with null variation of both constants. However, the results do not exclude the possibility that the 
fundamental constants have values different from their present ones but constant in the early universe. It is possible 
to obtain a linear relationship (between ^ and from the BBN and CMB contours: 

^ ^ ^ ^ = aBBN— + bBBN for BBN, (9) 

< V >o ao 

^ ^ ^ ^ = acMB— + bcMB for CMB, (10) 

< V >Q ao 

where qbbn = 0.181 ± 0.003, Bbbn = 0.0046 ± 0.0002, gcmb = -3.71^:5 and bcMB = -0.053tomT 

Figure [T^ show the 2D contours for — and .'^"V obtained from BBN and CMB data. In this case, the correlation 
coefficients are —0.82 for CMB and —0.91 for BBN. A phenomenological relationship between the variation of the 
fundamental constants a and mg can be obtained by adjusting a linear function. These two linear fits are different 
for both cases: 

7 ^ = CBBN ^dBBN for BBN, (11) 

(me)o "0 

7 f = ccMB + dcMB for CMB, (12) 

where cbbn = -1.229 ± 0.008, (Ibbn = -0.0234 ± 0.005, ccmb = and dcMB = -0.053^[!;j]^? (the time 

variation of < w > during CMB has the same effects than the variation of the electron mass) . 

It is important to point out that BBN degeneration suggests phenomenological relationships between the variations 
of both constants, while the CMB contours are not thin enough to assure any conclusion. 

Our results suggest that the model used by Ichikawa et al. where the variation of fundamental constants is 
driven by the time evolution of a dilaton field can be discarded, since these models predict me ~ a^/^. 
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FIG. 11: 2D contour levels for variation of a and < v > from BBN (solid line) and CMB (dotted line) data. 
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FIG. 12: 2D contour levels for variation of a and me from BBN (solid line) and CMB (dotted line) data. 

V. SUMMARY AND CONCLUSION 

In this work we have studied the joint time variation of the fine structure constant and the Higgs expectation value 
and the joint variation of a and me in the early Universe. We used the observational abundances of D, ^He and 
'^Li to put bounds on the joint variation of a and < t; > and on the joint variation of a and rUe during primordial 
nucleosynthesis. Wc used the three year WMAP data together with other CMB experiments and the 2dfGRS power 
spectrum to put bounds on the variation of a and < v > (or me) at the time of neutral hydrogen formation. 

Prom our analysis we arrive to the following conclusions: 

1. The consideration of different values of leads to similar constraints on the time variation of the fundamental 
constants. 

2. We obtain non null results for the joint variation of a and < v > at 6a in two cases: i) when tjb is allowed to 
vary and all abundances are included in the data set used to perform the fit and ii) when only ^He and ^Li are 
included in the data set used to perform the fit and r]B is fixed to the WMAP estimation. In the first case, the 
obtained value of ?7s is inconsistent with the WMAP estimation. 

3. We obtain non null results for the joint variation of a and mg at 6(7 in two cases: i) when rfg is allowed to 
vary and all abundances are included in the data set used to perform the fit and ii) when only ^He and ^Li are 
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included in the data set used to perform the fit and rjE is fixed to the WMAP estimation. In the first case, the 
obtained value of rjB is inconsistent with the WMAP estimation. 

4. We also obtain non null results for the joint variation of a and < v > and a and nie when all abundances are 
included in the data set and tjb is fixed to the WMAP estimation. The statistical significance of these results 
is too low to claim a variation of fundamental constants. 

5. Excluding ^Li abundance from the data set used to perform the fit, and keeping t]b fixed, we find results that 
are consistent with no variation of fundamental constants within la. 

6. The bounds obtained using data from CMB and 2dFGRS are consistent with null variation of a and < t; > (or 
nie) at recombination within la. 

7. We find phenomenological relationships for the variations of a and < v >, and for the variations of a and TOg, at 
the time of primordial nucleosynthesis and at the time of recombination. All the phenomenological relationships 
correspond to linear fits. 

8. From our phenomenological approach, it follows that the relationship between the variations of the two pairs of 
constants considered in this paper is different at the time of nucleosynthesis that at the time of neutral hydrogen 
formation. 

9. The dilaton model proposed by Ichikawa et al. [11] can be discarded. 
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APPENDIX A: PHYSICS AT BBN 



We discuss the dependencies on a, < u >, and TOg of the physical quantities involved in the calculation of the 
abundances of the light elements. We also argue how these quantities are modified in the Kawano code. 



1. Variation of the fine structure constant 



The variation of the fine structure constant affects several physical quantities relevant during BBN. These quantities 
are the cross sections, the Q-values of reaction rates, the light nuclei masses, and the neutron-proton mass difference 
(along with the neutrons and protons initial abundances and the n <-> p reaction rates). 

The cross sections were modified following refs. [4l|, [H, HI, Ull and replacing a by ao (^l + in the numerical 

code. The Q-values of reaction rates where modified following ref. [ill . 

To consider the effect of the variation of the fine structure constant upon the light nuclei masses, we adopted: 

Am,. „Aa , . , 

= P— , (Al) 



where P is a constant of the order of lO^'' (see ref. 4J] for details) and rrix is the mass of the nuclei x. These changes 
affect all of the reaction rates, their Q-values and their inverse coefficients. 

If the fine structure constant varies with time, the neutron- proton mass difference also changes. Following ref. [lOll |: 

^ -0.587^ . (A2) 
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This modifies the n <^ p and the neutrons and protons initial abundances. The n ^ p reaction rate is calculated by: 



An- 



K 



EePe 



{Ee + Am„p) 



K 



dE,^ 



EePe 



{Ee - Am„p) 



(A3) 



where if is a normalization constant proportional to G^p, E^ and Pe are the electron energy and momentum re- 
spectively, T~f and are the photon and neutrino temperature and is the ratio between the neutrino chemical 
potential and the neutrino temperature. In order to include the variation of Ammp we replace this quantity by 



( 1 I £Am^\ 



in Kawano code. The neutrons and protons initial abundances are calculated by: 

1 



I _^ gAm„p/T+C ' 
1 

1 _|_ g-Am„p/T-^ 



(A4) 
(A5) 



Variation of the electron mass 



If the electron mass can have a different value than the present one during primordial nucleosynthesis, the sum 
of the electron and positron energy densities, the sum of the electron and positron pressures and the difference of 
the electron and positron number densities must be modified in order to include this change. These quantities are 
calculated in Kawano code as: 



Pe- + Pe+ 
P. 



2 {mec^y 
^ {hcf 



^(-l)"+icosh (n0e)M(r 



+ 2 (m,c2) (^)»+i , 
— cosh {n(j)e) N(nz) 



[hcf 



TT 

Y 



n 3 



'^(-l)"+isinh (n0e)L(nz) 



(A6) 
(A7) 
(A8) 



where z = -^ffr^, 0e is the electron chemical potential and Liz), M[z) and A^(z) are combinations of the modified 

Bessel function Ki{z) [8llll02l |. The change in these quantities affects their derivatives and the expansion rate through 
the Friedmann equation: 



Stt 



A 



where G is the Newton constant, A is the cosmological constant and 

PT = P-r + Pe- + Pe+ + Piy+ Pb 



(A9) 



(AlO) 



The n p reaction rates (see Eq. (jA3|) ) and the weak decay rates of heavy nuclei are also modified if the electron 
mass varies with time. 

It is worth while mentioning that the most important changes in the primordial abundances (du e to a change in 
TTig) arrive from the change in the weak rates rather than from the change in the expansion rate 



3. Variation of the Higgs vacuum expectation value 



If the value of < w > during BBN is different than the present value, the electron mass, the Fermi constant, the 
neutron-proton mass difference and the deuterium binding energy take different values than the current ones. The 
electron mass is proportional to the Higgs vacuum expectation value, then 



Ame 
(me)o 



A < w > 

< V >o 



(All) 



16 



The Fermi constant is proportional to < v > [l03| |: this dependence affects the n ^ p reaction rates. The neutron- 
proton mass difference changes by [lOlj 

- — ^ = 1.587 , A12 

Aninp < V >o 

affecting n p reaction rates (see Eq. ljASp ) and the initial neutron and proton abundances (see Eq. (|A4p '). 
The deuterium binding energy must be corrected by 

^ = , (A13) 

where k is a model de pend ent constant. This constant can be found: i) using chiral perturbation theory, as was done 
by Beane and Savage |l04| : ii) using effective potentials to describe the nucleon-nucleon interaction. This correction 
affects the initial value of the deuterium abundance 

Yd = ^ , (A14) 

0.47110-iorg/^ 

where Tg is the temperature in units of lO^K, and ejj is in MeV. 



APPENDIX B: PHYSICS AT RECOMBINATION 



During recombination epoch, the ionization fraction, Xe — n^jn (where ng and n are the number density of 
free electrons and of neutral hydrogen, respectively), is determined by the balance between photoionization and 
recombination. 

In this paper, we solved the recombination equations using RECFAST ^99^, taking into account all of the depen- 
dencies on a and me [13, [H, [13, [s^ • To get a feeling of the depend encie s of the physical quantities relevant during 
recombination, we consider here the Peebles recombination scenario (l05| . The recombination equation is 



'dt \ n) 



^^|^a^_^^r^^_(s,-B.)/feT^ ^ (Bl) 



where 



(l + X(/3, + A2«,i.)ni,) ^ ' 

is the Peebles factor, which inhibits the recombination rate due to the presence of Lyman-a photons, riig is the 
number density of hydrogen atoms in the ground state, and -B„ is the binding energy of hydrogen in the nth principal 

quantum number. The redshift of the Lyman-a photons is if = with Aq. = and A2s,is is the rate of decay 

of the 2s excited state to the ground state via 2-photon emission, and scales as a^nie- Recombination directly to the 
ground state is strongly inhibited, so the case B recombination takes place. The case B recombination coefficient etc 
is proportional to a^irie The photoionization coefficient depends on ckc, but it also has an additional dependence 
on me, 

P. = aJ^^)"\--^'-- . (B3) 



The most important effects of changes in a and me during recombination are due to their influence upon Thomson 
scattering cross section ut — ^^2^2 o? , and the binding energy of hydrogen Bi = ^ a^meC^. 
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